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We consider two perturbative schemes to calculate excitation energies, each employing the Kohn- 
Sham Hamiltonian as the unperturbed system. Using accurate exchange-correlation potentials gen- 
erated from essentially exact densities and their exchange components determined by a recently 
proposed method, we evaluate energy differences between the ground state and excited states in 
first-order perturbation theory for the Helium, ionized Lithium and Beryllium atoms. It was re- 
cently observed that the zeroth-order excitations energies, simply given by the difference of the 
, Kohn-Sham eigenvalues, almost always lie between the singlet and triplet experimental excitations 

energies, corrected for relativistic and finite nuclear mass effects. The first-order corrections pro- 
vide about a factor of two improvement in one of the perturbative schemes but not in the other. 
^ ■ The excitation energies within perturbation theory are compared to the excitations obtained within 

ASCF and time-dependent density functional theory. We also calculate the excitation energies in 
perturbation theory using approximate functionals such as the local density approximation and the 
optimized effective potential method with and without the Colle-Salvetti correlation contribution. 
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I. INTRODUCTION 



e 

73 ; 

£^ , Density functional theory is usually employed to study ground states. However, it is clear that density functional 
theory can in principle be used to calculate excitation energies since, according to the Hohenberg-Kohn theorem [[l], 
the density determines the external potential (aside from an additive constant) and, consequently, the entire spectrum. 

^ . In principle, one could construct functionals of the ground state densities that yield excitation energies. 

•i-H . 

There are several proposals in the literature to deal with excited states. Gunnarson and Lundquist |^j pointed out that 
the original Hohenberg-Kohn theorems can be generalized to apply to the lowest energy state of any given symmetry. 
Von Barth || extended density functional theory to states of mixed symmetry following the observation of Ziegler et 
al. [Q that the Xcv method is designed to describe states that have single determinant wave functions in the absence 
of interaction. The practical recipe consists of calculating, within the local density approximation, the energies of 
states of mixed symmetry represented by a single determinant and linearly combining them to construct the energies 
of the states of pure symmetry Further developments and applications of this procedure, known as ASCF, can 

be found in Refs. Another generalization of density functional theory to treat excited states is the ensemble 

density functional theory |l0|-[l^1 . Recently, a different and quite promising method to calculate excited states based 
on time-dependent density functional theory has been developed p^3-n7[| . 



An alternative approach [ fl8| is to calculate excitation energies perturbatively using the non-interacting Kohn-Sham 
Hamiltonian as the zeroth-order Hamiltonian. In the Kohn-Sham formulation of density functional theory Jl9| , the 
interacting electron system is replaced by a system of non-interacting electrons in an effective local potential that 
yields the true ground state density. In the perturbative approach, the zeroth-order excitation energies are simply 
the differences of Kohn-Sham eigenvalues. Remarkably, it has been recently shown ]29] that, for atomic systems, the 
excitation energies obtained as differences of Kohn-Sham eigenvalues lie between the triplet and singlet experimental 
excitation energies, corrected for relativistic and finite nuclear mass effects, with the exception of the Is — > 3d 
transition of ionized Lithium. This finding is surprising since, except for the highest occupied orbital, there is no 
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strict relationship between the Kohn-Sham eigenvalues and excitation energies. 



In an attempt to improve on these results and determine the singlet-triplet splitting, we calculate excitation energies 
to first-order in perturbation theory. We employ two variants of perturbation theory. Interestingly we find that 
the first-order corrections from the more obvious standard perturbation theory make the excitation energies worse on 
average while those from the coupling- constant perturbation theory of Gorling and Levy | pl| yield about a factor of two 
improvement. To perform these calculations, we employ excellent approximations to the true Kohn-Sham exchange- 
correlation potentials determined by generating an accurate density for the system of interest and then computing 
an exchange-correlation potential that yields the desired density as the ground state solution for the fictitious non- 
interacting system |22|-|24| . Evaluating the first-order corrections also requires the knowledge of the correct exchange 
component of the exchange-correlation potential as determined in our previous work | p5[ . Standard perturbation 
theory has also been recently applied to calculate ground state energies and the perturbation series based on various 
approximate exchange-correlation functionals were found to be divergent for the systems tested [E6j . 

In Sec. |, we briefly introduce density functional theory and its Kohn-Sham formulation. In Sec. |l| we present the 
coupling-constant scheme and the formulae used to determine the excitation energies to first- order in the coupling- 
constant. In Sec. Ill, we give the results and a comparison with other methods to evaluate excitation energies. The 
perturbation theory scheme is also implemented using approximate exchange-correlation functionals. In Appendix |a|, 
we derive the formulae for the first-order corrections to the excitation energies. 



II. THEORETICAL BACKGROUND 



In the Kohn-Sham formulation of density functional theory 1X91] , the ground state density is written in terms of 
single-particle orbitals obeying the equations in atomic units (fi = e = m = 1): 

j-^V 2 + Wext(r) + vn(r) + v xc {[p] ; r)| ^ = e^i, (1) 

where 

N 

p(r) = ^|^(r)| 2 . (2) 

i=l 

The electronic density is constructed by summing over the JV lowest energy orbitals where iVis the number of electrons. 
v ex t( r ) is the external potential and vn(r) the Hartree potential. The exchange-correlation potential v xc ([p] ;r) is the 
functional derivative of the exchange-correlation energy E xc [p] that enters in the expression for the total energy of 
the system: 

E = -\j2j ^V 2 ^ dr + J p (r) v cxt (r) dr + \ J J P ^ )p{ ^ dr dr' + E xc [p] . (3) 

The exchange-correlation functional is written as the sum of two separate contributions for exchange and correlation, 

E xc [p] = E x [p] + E c [p] , (4) 

yielding a corresponding splitting in the exchange-correlation potential, v xc ([p);r) = v^([p\; r) + v c ([/?]; r). The defini- 
tion of the exchange energy is in terms of the non- interacting wave function <E>0i the Slater determinant constructed 
from the Kohn-Sham orbitals, as 

E x [p] = (# | V cc \*o)~\Jj ^^dr dr', (5) 
where V oc is the electron-electron interaction. 
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III. PERTURBATION THEORY 



In many-body perturbation theory the Hamiltonian TL is written as the sum of a zeroth-order Hamiltonian TLo and a 
perturbing Hamiltonian TL\ = TL — TLq. 

In the commonly used Moller-Plesset perturbation theory [^7|, the zeroth-order Hamiltonian is the Fock operator: 

H% F = T + V cxt + Vh + K HF , (6) 

while in density functional perturbation theory 

H = T + V eKt + V B + V m) (7) 

where T is the kinetic operator and V^ HF is the non-local Hartree-Fock exchange operator. Here and elsewhere, we em- 
ploy the notation that potentials are capitalized when summed over all the electrons. The perturbation corresponding 
to the zeroth-order (Eq. 0) is clearly given by 

Hi - V„ — Vh — V xc . (8) 

The other variant of perturbation theory we consider is due to Gorling and Levy pl| ] and is based on the adiabatic 
connection Hamiltonian which links the Kohn-Sham Hamiltonian with the fully interacting Hamiltonian keeping the 
ground state density constant independent of a: 

H a = T+ K xt + V a + V KC + a (V cc -V x - Vh) - V c a . (9) 

v™ is second order in a PH and equals the correlation potential when a = 1. The zeroth-order Hamiltonian is again 
the Kohn-Sham Hamiltonian and the perturbing Hamiltonian contains a term 

H[ 1] = a (V ee -V s - ty , (10) 

which is linear in the perturbation parameter a and a component — v" which contains second and higher order 
contributions pjfl . To any order in perturbation theory, the density equals the true ground state density to that order 
in the coupling-constant [|TJ. 

Let us consider the excitation of an electron from the space-spin orbital <j) k to <f> v . In both perturbative schemes, 
the zeroth-order energy difference equals the difference in the expectation value of the Kohn-Sham Hamiltonian on 
the eigenstates $ and $fc^, obtained from $ by exciting an electron from orbital k to orbital and is simply the 
difference of the Kohn-Sham eigenvalues: 

AE (°) = e v - e fc . (11) 

The correction to first-order in the coupling-constant (Eq. ^) is given by 

AE« = \V CC - Vh - V x \ - ($ \V CC - V H - V x \ $> . (12) 

This expression simplifies for the Helium, ionized Lithium and Beryllium atoms, that we use as examples in this work, 
all having ground states of 1 S symmetry. The excited states are either singlets or triplets. Let $0 be the Kohn-Sham 
ground state and the excited state obtained by exciting an electron from a (doubly occupied) orbital k to orbital 
v. In this case, as shown in Ref. |Q and Appendix the first-order energy difference between the ground state and 
the state <& kv depends on whether $> ku is a triplet or a singlet and is given respectively by 

AE« (T, k -> v) = (fa \v^ F - v x \ fa) - (fa \v* F - v x \ fa) - (V»„V* \g\ Mk) , (13) 

and 

AE« (S, k v) = (fa \v^ F - v x \ fa) - (fa \vT - v x \ fa) - {M k \g\ ^ k ) + 2 {^fa \g\ fa$ v ) , (14) 

where g= l/|r — r'| and if> is just the spatial component of the orbital </>. is the Hartree-Fock exchange potential, 
but constructed from the Kohn-Sham orbitals occupied in ground state ^o- 

To evaluate the correction within the standard perturbation theory (Eq. we simply replace v x with v xc in the 
above formulae. 
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IV. RESULTS AND DISCUSSION 



In order to evaluate the excitation energies to first-order in the coupling constant (Eqs. [TT], [ll| and |lj) and in the 
more standard perturbation scheme, we need to determine an accurate exchange-correlation potential, Kohn-Sham 
orbitals and eigenvalues and the exchange component of the exchange-correlation potential for the system of interest. 

As mentioned in Sec. [l| the accurate exchange-correlation potential used in Eq. |l| for the determination of the Kohn- 
Sham orbitals and eigenvalues is obtained starting from an accurate density and searching for a potential that yields 
the desired density as a solution of the Kohn-Sham equations (Eq. |l| and ||) . In the special case of the singlet ground 
state of a two-electron system [p8pE(| , the exchange-correlation potential can be simply obtained from Eq. |l]. For 
systems with more than two electrons, v xc can be determined by expanding it in a complete set of basis functions and 
varying the expansion coefficients such that Eqs. [I] and |^ yield the accurate density (22 - 24 1 . The exchange component 



of the exchange-correlation potential is simply given by the condition that it cancels the self-interaction term in the 
Hartree potential in the case of two electrons in the singlet ground state. For a many-electron system, the exchange 
potential is obtained as discussed in Ref. p5| . 

In Fig. and in Tables @, [n] and H|, we show the excitation energies to zeroth-order and to first-order in the two 
perturbative schemes for Helium, ionized Lithium and Beryllium respectively. 

The zeroth-order excitation energies are simply given by the difference in the Kohn-Sham eigenvalues and, as already 
discussed in Ref. pc|| , lie between the experimental triplet and singlet excitation energies to an accuracy of four 
significant digits, provided that relativistic and finite nuclear mass effects are taken into account. The only exception 
is the Is — * 3d transition for ionized Lithium. 

The first-order excitation energies obtained from standard perturbation theory are on average less accurate than 
the zeroth-order energies. In contrast, the excitation energies to first-order in the coupling-constant are on average 
better than the zeroth-order energies but the improvement is not consistently present for all the excitations. With 
the exception of the 2s to 3d, Ap and Ad excitations for Beryllium, the singlet-triplet splittings, obtained from cither 
first-order scheme, are too large. Note that for the two-electron systems, as shown in Appendix |A|, Eqs. [l3] and |lj 
yield a symmetrical splitting around the zeroth-order excitation energy. 

From Fig. |l|, it is clear that the theoretical predictions for the higher excitation energies converge to the experimental 
excitation energies both in zeroth-order and in first-order coupling constant perturbation theory. On the other hand, 
the excitation energies from first-order standard perturbation theory converge to a too large value for Helium and 
ionized Lithium and a too low value for Beryllium. It is apparent that the zeroth-order energies must converge to the 
experimental energies since the highest occupied Kohn-Sham eigenvalue is the negative of the ionization energy [j30| . 
In Fig. |l|, the ionization limit for each system has been indicated by an arrow. There is an explanation of why 
the higher excitation energies from first-order coupling constant perturbation theory converge to the experimental 
excitation energies [ fn| . Coupling constant perturbation theory is based on the coupling constant Hamiltonian in 
which the external potential is varied to keep the charge density constant for all values of the coupling constant a. 
Since the ionization energy is related to the asymptotic decay of the charge density, it is independent of the coupling 
constant. Now, let us consider a perturbation expansion of the ionization energy and charge density in powers of the 
coupling constant. Since the ionization energy and density are independent of the coupling constant, each term in the 
series, except the zeroth-order term, must be zero. If we perform perturbation theory to a finite order, the ionization 
energy is still strictly constant at each order: since the ionization energy is obtained as an energy difference, each order 
of perturbation theory contributes only to that same order. We note in passing that this is not true of the density. 
At a given order of perturbation theory, the density is correct only to that order since the n-th order wavefunction 
gives an (n + l)-th order contribution to the charge density that is canceled only by going to the (n + l)-th order of 
perturbation theory. Finally, we note that the reason why the first-order energies in the standard perturbation theory 
are too high for Helium and ionized Lithium and too low for Beryllium is that the correlation potential, v c , is mostly 



positive for the former and entirely negative for the latter |2J|. As explained in Sec. [II, the first-order excitation 
energies in standard perturbation theory differ from the energies to first-order in the coupling constant in having v xc 
instead of only v c in Eqs. [ll] and |TJ. 

We observe from Fig. [| that the zeroth-order and first-order in coupling constant results are accurate not only close to 
the ionization limit but over most of the energy range. This calls for some further explanation beyond that provided 
above. In Ref. the accuracy of the zeroth-order excitation energies is attributed to the additional fact that the 
Kohn-Sham orbitals and the quasi-particle amplitudes satisfy the same equation to order 1/r 4 . 

We wish to emphasize that the accuracy of the first-order excitation energies, obtained as differences of first-order 
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total energies of the ground and excited states, does not extend to the individual total energies of the ground and 
excited states. Prom Eq. |[ we obtain for the ground state energy to first-order in the coupling constant: 

E« = ($o|ft (0) +ft (1) |$o> = (4o|T + V r ox t+Veo + V r c|*o>. (15) 

Now, ($o\T + V cx t + V ee |$o) yields an energy that is a little worse than the energy from either the optimized 
effective potential (OEP) method J3^] or the Hartree-Fock method since the Hartree-Fock wavefunction is the best 
single-determinant wavefunction, while ($o| V c \ &o) can be either negative or positive. For example, it is negative for 
Beryllium but positive for Helium, ionized Lithium and Neon. Consequently, at least for some systems, the ground 
state energies obtained from the first-order formula are worse than those from Hartree-Fock and the typical errors in 
Hartree-Fock total energies is considerably larger than the errors we obtain for the perturbative excitation energies. 



We summarize the results in Table IV which shows the root mean square errors in the excitation energies for our 
three test systems in zeroth-order and first-order perturbation theory. For the perturbation in the coupling-constant, 
the errors are reduced by about a factor of two for Helium and Beryllium and a factor of three for ionized Lithium. 
For the standard perturbation theory, the errors are about a factor of two worse than in zeroth-order for Helium and 
Beryllium while a slight improvement is obtained for ionized Lithium. 

In addition to the excitations from the highest occupied level, that we have considered so far, we also studied the 
core-level excitation of Beryllium from the Is orbital to the 2p orbital, Is 2 2s 2 1 S — ► Is 1 2s 2 2P 1 X P. The experimental 
energy is 4.243 Hartree p3| , while the zeroth-order perturbative energy is 4.017 Hartree and the first-order perturbative 
energies for the 3 P and H^ 1 states are 4.440 and 4.460 Hartree respectively. Clearly, neither the zeroth- order nor the 
first-order energies agree with experiment as closely as do the energies of excitations from the highest occupied level. 
This is related to the fact that the eigenvalues of the occupied orbitals, other than the highest, are not strictly related 
to the corresponding ionization energies. 

Just as Moller-Plesset perturbation theory is commonly used in quantum chemistry calculations, density functional 
functional perturbation theory could become a feasible technique to access excited states. However, the results 
presented above were obtained using the true Kohn-Sham potentials whereas practical implementations of the method 
would need to use potentials obtained from approximate exchange-correlation functionals. In Fig |^, we compare the 
zeroth-order excitation energies obtained from the local density approximation (LDA) and the optimized effective 
potential method with those obtained from the true Kohn-Sham potential. The LDA does not bind any of the 
unoccupied orbitals in Helium and binds only the lowest two unoccupied orbitals in Beryllium. Hence, the sparsity of 
LDA data in Fig ||. Since the highest occupied level in LDA is always too shallow, the zeroth-order LDA excitation 
energies are on average too small. However, note that the lowest excitation energies of Beryllium from all three 
potentials are very close. The reason is that this excitation is within the same principal quantum number shell 
(2s — > 2p) and, consequently, both orbitals involved in the excitation explore the same, rather limited, spatial region. 
The OEP excitation energies are too high for Helium and ionized Lithium and too low for Beryllium, which of course 
is due to the correlation potential being predominantly positive for the former and negative for the latter (24) . 

In Fig. |^, we present the same comparisons for the excitation energies obtained within first-order coupling-constant 
perturbation theory. The first-order OEP and first-order exact Kohn-Sham energies split in a roughly symmetrical 
way about the corresponding zeroth-order values while the first-order LDA energies have a large positive shift relative 
to the LDA zeroth-order. The reason is the following. The second term in Eqs. O and [l4|, — (0^ |t>^ F — uJ 4>k), is 
strictly zero for both the true Kohn-Sham potential |3(| and for the OEP potential pTI] when <pk is the highest occupied 
level. On the other hand, this identity does not hold for the LDA exchange potential which is in fact considerably 
shallower than the true one, resulting in a large positive contribution from this term. Note that all three quantities, 
<f>k, v^ F and v K , depend on whether they are being calculated for the true Kohn-Sham theory, OEP or LDA but the 
main difference comes from the change in u x . 

We do not show GGA results in Figs. || and || but we find that the existing GGA's do not yield a large improvement 
upon the LDA results. This is to be expected since the GGA potentials also fail to have the correct asymptotic decay 
at large distances and are too shallow for much of the region where the charge density is appreciable. 

Recently, there has been considerable interest in using time-dependent density functional theory (TDFT) to calculate 
excitation energies. In performing these calculations, a choice has to be made for the static potential and for the 
dynamic response. In Fig. |i[ we compare the TDFT excitation energies of Helium, obtained by Petersilka, Gossman 
and Gross JlTj using our exact Kohn-Sham potential and the dynamic response from OEP, to our first-order coupling- 
constant perturbation theory results. The agreement is very close and on the scale of the plot the difference between 
the two theories can only be detected for the lowest two excitations. For two electron systems, it is indeed possible 
to show that the TDFT excitation energies obtained within the single-pole approximation, using the exact Kohn- 



5 



Sham potential and the time-dependent OEP kernel, are strictly identical to the excitation energies derived in first- 
order perturbation theory in the coupling constant (Eqs. [l3| and |l4|). The fact that, in this case, the single-pole 
approximation gives the leading correction to the Kohn-Sham eigenvalue differences |l7j explains the close agreement 
found in Fig. ^. 

As mentioned in the Introduction, an older method for calculating excitation energies is the ASCF method 
Relative to perturbation theory, ASCF has the advantage that, when using approximate functionals, it can be applied 
even when the corresponding approximate exchange-correlation potential does not bind any unoccupied orbital. In 
Fig. [| we compare the excitation energies of Beryllium from ASCF, using the LDA potential. The results obtained 
for Beryllium are somewhat more accurate than those from perturbation theory using the LDA or OEP potentials 
but of course not as good as the perturbation theory results using the true Kohn-Sham potential. Also for ionized 
Lithium, we find that, although less accurate than the perturbation results with the exact Kohn-Sham solutions, the 
ASCF energies are more accurate than the LDA perturbative energies and of comparable accuracy to the perturbative 
OEP results. 

One can ask if it is possible to improve upon the perturbation theory results obtained with the exact exchange (OEP) 
functional by adding an approximate correlation functional. Unfortunately, the correlation potentials obtained from 
the LDA and various GGA correlation functionals bear no resemblance to the true correlation potential |Q . Recently, 
there has been considerable interest pql in orbital-dependent correlation functionals, the most popular one being the 
Colle-Salvetti functional . In Figs - ]!] and ||, we also show the correlation energies obtained from adding the Colle- 
Salvetti correlation functional to the true exchange (OEP) functional. The energies are improved for Beryllium but 
are worse for Helium and ionized Lithium since the Colle-Salvetti correlation potential is mostly negative for all three 
systems whereas the true correlation potential is negative for Beryllium but mostly positive for Helium and ionized 
Lithium. 

Clearly, a present obstacle to the success of perturbation theory in the coupling constant seems to be the lack of 
an accurate approximation for the exchange-correlation functional. Since approximate functionals such as the local 
density or generalized gradient approximations bind few (if any) unoccupied orbitals, combining an exact treatment 
of exchange with a local orbital-dependent correlation functional, better than the Colle-Salvetti functional, appears 
to be the most promising possibility in this direction. 

TABLE I. Excitation energies of He in Hartree atomic units. The theoretical energies of Drake and the perturbative 
results are for infinite nuclear mass and neglect relativity. 0(0) is the zeroth-order perturbative result. 0(1) and O a (l) are 



the first-order results 


in the standard and 


in coupling- const ant perturbation 


schemes respectively. 






Transition 


Final State 


Experiment 


Drake 


O(0) 


0(1) 


Oc(l) 


Is -> 2s 


2 3 S 
2 1 S 


0.7283 
0.7576 


0.7285 
0.7578 


0.7460 


0.7374 
0.7829 


0.7232 
0.7687 


Is -> 2p 


l 3 Po 

I'Po 


0.7704 
0.7797 


0.7706 
0.7799 


0.7772 


0.7833 
0.7989 


0.7693 
0.7850 


Is — > 3s 


3 3 S 
3 X S 


0.8349 
0.8423 


0.8350 
0.8425 


0.8392 


0.8483 
0.8594 


0.8337 
0.8448 


Is — > 3p 


2 3 P a 
2 1 P a 


0.8455 
0.8484 


0.8456 
0.8486 


0.8476 


0.8598 
0.8645 


0.8453 
0.8500 


Is -> 3d 


1 3 D 
I'D 


0.8479 
0.8479 


0.8481 
0.8481 


0.8481 


0.8629 
0.8630 


0.8481 
0.8482 


Is 4s 


4 3 5" 
4 X S 


0.8670 
0.8700 


0.8672 
0.8701 


0.8688 


0.8814 
0.8857 


0.8667 
0.8710 



G 



TABLE II. Excitation energies of Li + in Hartree atomic units. The theoretical energies of Drake pq] and the perturbative 



results are for infinite nuclear mass and neglect relativity. 0(0) is the zeroth-order perturbative result. 0(1) and O a (l) are 
the first-order results in the standard and in coupling- const ant perturbation schemes respectively. 



Transition 


Final State 


Experiment 


Drake 


0(0) 


0(1) 


Q (1) 


Is -> 2s 


2 3 S 
2 1 S 


2.1690 
2.2330 


2.1692 
2.2390 


2.2082 


2.1748 
2.2655 


2.1628 
2.2535 


Is -> 2p 


l 3 Po 


2.2521 
2.2865 


2.2522 
2.2866 


2.2733 


2.2612 
2.3079 


2.2499 
2.2966 


Is — > 3s 


3 3 S 
3 l S 


2.5277 
2.5460 


2.5278 
2.5462 


2.5375 


2.5387 
2.5615 


2.5261 
2.5489 


Is — > 3p 


2 3 P a 
2 1 P 


2.5493 
2.5595 


2.5495 
2.5597 


2.5554 


2.5611 
2.5745 


2.5487 
2.5621 


Is -> 3d 


1 3 D 
1 1 D 


2.5572 
2.5574 


2.5574 
2.5575 


2.5576 


2.5702 
2.5706 


2.5574 
2.5578 


Is — > 4s 


4 3 5" 
4 X S 


2.6426 
2.6499 


2.6428 
2.6501 


2.6466 


2.6548 
2.6638 


2.6421 
2.6511 
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TABLE III. Excitation energies of Be in Hartree atomic units. The perturbative results are for infinite nuclear mass and 



neglect relativity. O(0) is the zeroth-order perturbative result. 0(1) and O a {l) are the first-order results in the standard and 
in coupling-constant perturbation schemes respectively. 



Transition 


Final State 




Experiment 


O(0) 


0(1) 


0a(l) 


2s -> 2p 


l 3 Po 
1 1 P 




0.1002 
0.1939 


0.1327 


0.0621 
0.1981 


0.0629 
0.1989 


2s -> 3s 


2 3 S* 
2 J 5 




0.2373 
0.2491 


0.2444 


0.2006 
0.2231 


0.2331 
0.2556 


2s -> 3p 


2 3 P D 
2 1 P 




0.2679 
0.2742 


0.2694 


0.2311 
0.2412 


0.2640 
0.2741 


2s -> 3d 


1 3 D 
1 x 73 




0.2827 
0.2936 


0.2833 


0.2460 
0.2498 


0.2814 
0.2852 


2s -> 4s 


3 3 S 
3 1 S 




0.2939 
0.2973 


0.2959 


0.2557 
0.2619 


0.2928 
0.2990 


2s -> Ap 


3 3 P D 
3 1 P a 




0.3005 
0.3063 


0.3046 


0.2660 
0.2691 


0.3029 
0.3061 


2s -> 4d 


2 3 D 
2 ^ 




0.3096 
0.3134 


0.3098 


0.2713 
0.2730 


0.3089 
0.3106 


2s -> 5s 


4 3 S 

/I 1 c 




0.3144 
0.3156 


0.3153 


0.2758 
0.2784 


0.3139 
0.3166 


TABLE IV. 
theory. 


Root mean square errors 


of the excitation energies 


of He, Li + and Be in zeroth-order and first-order perturbation 


System 




0(0) 




0(1) 




0a(l) 


He 
Li+ 
Be 




0.0067 
0.0167 
0.0178 




0.0159 
0.0148 
0.0358 




0.0039 
0.0056 
0.0099 
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S 0.28 



.2 0.20 



j3 0.08 



Beryllium: perturbation 


theory 


o 0(0) 
+ ©(1) 
* 6>a0) 


+ 

/ + 
+ 


/ 









0.04 0.08 0.12 0.16 0.20 0.24 0.28 0.32 0.36 
Experimental excitation energy (Hartree) 



FIG. 1. Theoretical versus experimental excitation energies for He, Li+ and Be (see Tables O(0) is the zeroth-order 

perturbative result. 0(1) and O a (l) are the first-order results in the standard and in coupling-constant perturbation schemes 
respectively. The ionization limits are indicated by arrows. 
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Experimental excitation energy (Hartree) 




0.08 0.12 0.16 0.20 0.24 0.28 0.32 

Experimental excitation energy (Hartree) 

FIG. 2. Zcroth order perturbation theory excitation energies for He, Li + and Be. The theoretical excitation energies obtained 
from the exact Kohn-Sham, LDA, OEP and OEP-Colle-Salvetti eigenvalues are plotted versus the experimental values. The 
singlet and triplet energies from the same configuration are joined by lines. 
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FIG. 3. Excitation energies for He, Li + and Be to first-order perturbation theory in the coupling constant. The theoretical 
excitation energies obtained from the exact Kohn-Sham solutions, LDA, OEP and OEP-Colle-Salvetti are plotted versus the 
experimental values. 
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Helium: perturbation theory 



W 0.85 



a 

.2 0.80 



o TDFT-OEP 



0.70 ^ 1 ' 1 1 

0.70 0.75 0.80 0.85 0.90 

Experimental excitation energy (Hartree) 



FIG. 4. Excitation energies for He within first-order perturbation theory in the coupling constant and time-dependent density 
functional theory. The time-dependent DFT excitation energies are obtained jl7| using the exact Kohn-Sham potential as static 
potential and OEP for the dynamical response. The theoretical excitation energies are plotted versus the experimental values. 




FIG. 5. Excitation energies for Be within first-order perturbation theory in the coupling constant and ASCF-LDA. The 
theoretical excitation energies are plotted versus the experimental values. 
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APPENDIX A: FIRST-ORDER CORRECTION 



We obtain here the expressions for the first-order energy differences given in Eqs. |l3] and We consider the ground 
state Kohn-Sham determinant $o and an excited determinant $fc„ obtained by exciting an electron from a doubly 
occupied orbital k to an empty orbital v. 

First, we rewrite the contributions to the first-order energy difference (Eq. |l2|) in terms of the Kohn-Sham orbitals: 

(§kv |«H + «x| ®kv) ~ ($0 |«H + Wx| $o) = {<j>u |«H + V x \ 4>v) - (0fc |«H + «x| 4>k) , (Al) 

and 

^ occ 

($\V ee \$) = -Y J ((i3\9\i3)-(ij\9\3i)), 

id 

where (ij \g\rs) = (0i0j \g\ <j> r <t> s ) — J J drdr'0*(r)0*(r')g0 r (r)0 s (r') and the sum is over the orbitals occupied in 
the determinant $. Since the Hartree potential i>h is constructed from the orbitals occupied in the ground state <£>0j 
we obtain 

occ 

\V ee \ $ fcv ) - ($0 \V ee \ * ) = £ ( - (kj \g\ kj) + (kj \g\ jk) + (uj \g\ vj) ~ (vj \g\ jv) ) - (vk \g\ vk) + (vk \g\ kv) 

3 

= — (k |<u H + «* F | k) + (v\v K + Vx F \ v) - (vk\g\vk) + (vk\g\kv) (A2) 

where the sums are over the Kohn-Sham orbitals o ccupied in <I>o and v^ F is the Hartree-Fock exchange potential 
constructed from these orbitals. Combining Eqs. [l2, Al , A2 , we obtain 









\€ F -v*\ 




0fc \Vx F - Vx\ 



(A3) 



Our derivation of Eq. A3 is for a single-determinant &kv If the final state is a triplet, its Kohn-Sham wave function 
can be written in terms of a single determinant provided that we consider either the S z = 1 or the S z = — 1 components 
of the triplet. Then, the third term is zero since we are considering the excitation from 0fc = ipkXl t° ^ = VvXt 
or from 0& = ipkXl ^° — ^vXi- If ^ku is a singlet, it is necessary to write as a sum of two determinants. 
There are two possible excitations, V'fcXT — > VVXT an d "0fcXi ~* 0VXJ.> an d the resulting wave function is equal to 
®kv = ($i - $2) /V2. Consequently, 

($k„ \V CC \ - \ [(^llKel^l) + (* 2 |Ke|*2) " (*l|Vee|4 2 ) - ($2 |K=c |*l)] • 



The first and second contribution are equal to each other and equivalent to the above calculation A3. The last two 
contributions are equal and, if we define 0i — ipkXli 02 = ipvXh 4>3 = V'feXT an d 4>i — VvXi, gi ye 

- i [(*i|V co |$ 2 ) + ($a|Vee|*l>] = - (01^ 2 |5l^3^4) + (0102 | 5 |0403> 

The first term is zero due to spin orthogonality while the second term can be rewritten as (ipkipv \g\ ?/v0fc)- Combining 
the above results, we obtain Eqs. [t3] and |lj. 

In Sec. [II, we mentioned that, for two-electrons in a singlet state, formulae (Eq. [l^ and |l4|) yield a symmetrical 
splitting around the zeroth-order excitation energy. The ground state Kohn-Sham determinant of two electrons in a 
singlet state is constructed from the two orbitals <f>i — V'oXt anc l 4>2 = 0oX|- Here, we consider the excitation of an 
electron from either 0x or 2 to the orbital 0„. 

Since, for this system, the exchange potential satisfies v x = —vn/2 = —{ipo\g\tl>o}, we obtain 

(0i| V X |0r) = (02 | Vx |02> = -(^oMalMo) = (01 1 Vx F I0l) = (02 I |0 2 ), 

so that (0x| v* F _ Wx = (0 2 | {)HF _ Ux |0 2 ) = 0) and 
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{<p v \v x - v x |<M = -{i^u^o\g\^o^ v ) + (VvV'ols'lVvV'o)- 

Consequently, by combining the above equations and the formulae (Eq. [l3]and [l4|) , we obtain a symmetrical splitting 
around the zeroth-order excitation energy: 

AE« (T,0 -f v) = - (^o|<?|V>oVv), 

AE« (S, -> i/) = ftMtoM^lM- (A4) 
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